ϵ-expansion of multivariable hypergeometric functions appearing in Feynman integral calculus

We present a new methodology, suitable for implementation on computer, to perform the ϵ-expansion of hypergeometric functions with linear ϵ dependent Pochhammer parameters in any number variables. Our approach allows one Taylor as well Laurent series expansion multivariable functions. Each coefficients is expressed combination same domain convergence that original function. illustrative examples one, two and three variables which are typical Feynman integral calculus.